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We study synchronization properties of general uncoupled limit-cycle oscillators driven by com- 
mon and independent Gaussian white noises. Using phase reduction and averaging methods, we 
analytically derive the stationary distribution of the phase difference between oscillators for weak 
noise intensity. We demonstrate that in addition to synchronization, clustering, or more generally 
coherence, always result from arbitrary initial conditions, irrespective of the details of the oscillators. 

Noise- induced synchronization is widely observed in various experimental systems ranging from neurons to lasers [l[ . 
From the theoretical standpoint, after severalpioneering studies [2|, significant progress has been made by utilizing the 
phase reduction method for limit cycles [3, 0, [j| . These works generally proved that when the limit-cycle oscillators are 
driven by a sufficiently weak common additive noise, the Lyapunov exponent of the synchronized state averaged over 
the noise distribution always becomes negative, namely, the synchronized state is statistically stabilized. However, 
these works are still incomplete as the Lyapunov exponent only characterizes local stability and do not describe global 
behavior of the oscillators. Also, effects of multiplicative common noises and non-vanishing additional independent 
noises remain unclarified. In this letter, we analyze this phenomenon in more detail from an alternative perspective by 
adopting phase reduction and averaging methods to many-body stochastic dynamical equations describing a general 
class of limit-cycle oscillators driven by common and independent noises, which yields global characterization of their 
synchronization properties. 

We consider the following Langevin equations describing an ensemble of N uncoupled identical oscillators driven 
by common and independent noises: 

X^ a \t) = F(X (Q >) + \/zX?(X (a) )£(t) + ^feH{X^)r]^{t), (1) 

for a = 1,---,A, where G R M represents the state of the a-th oscillator at time t, G R M its 

individual dynamics, £(t) G R M the external noise common to all oscillators, and ri^ a >{t) £ R M the external noise 
added independently to each oscillator. and r}( a > (i) are assumed to be independent, identically distributed 

zero-mean Gaussian white noises of unit intensity and correlation functions given by (£i(t)£j(s)) = Si.j5(t — s), 

(vl a) (tWf\s)} = 6 a ,fi5ij5(t-8), and (^(t)^ a) (s)) = (the subscript i or j denotes the vector component). The 
parameters D and e control their intensities. The R MxM matrices G(X^) and H^X^) represent the coupling 
of the oscillator to both types of the noises, which are assumed to be smooth functions of X^ a \ We interpret these 
Langevin equations in the Stratonovich sense, namely, we consider the white noise as the limit of colored noise with 
vanishingly small correlation time. 

We assume that (i) each oscillator obeys the same dynamics, with a single stable limit cycle Xnft) in its phase 
space [lH , and that (ii) noises of both types are sufficiently weak, so that phase reduction 0, 0, Q of the above 
Langevin equations is possible [l3| . Specifically, we describe the dynamics of each oscillator using only a constantly- 
increasing phase variable <p(i) = 4>(X(t)) 6 [— 7r,7r], defined along its limit cycle and also on its phase space except at 
phaseless sets. Applying the standard phase reduction method to Eq. ((T|) [||, we obtain (by virtue of the Stratonovich 
interpretation) the following approximate Langevin equations for the phase variables cf> = (t^ 1 -*, • • • , ^ N ')' 

0<°>(t) = u + VDZ(4>^) • G(0 (Q) )£(t) + V~eZ((f> (a) ) ■ H(4>^) V ^(t). (2) 

Here, u> is the natural frequency of the oscillators, Z((p a ') — V x4>" a ^ lx=x (>(°>) R M ^ s the phase sensitivity 
function of the individual oscillator that quantifies the phase response of each oscillator to weak perturbations 0], 
G{(f>^) = G(X (</>( tt ))), and H{<j>^) = H(X (<f>^)). We normalize Z{<j>) such that Z(<p) ■ F{X (<f>)) = u> holds 
constantly. Z (</)), G((f>), and H(4>) are smooth periodic functions of cj>. 

The Stratonovich Langevin equations ^ are converted to equivalent Ito stochastic differential equations [l(| of 
the form d(f>( a \t) — [<p)dt + d^ a ^ [cj>, t) , where {(^ (cf>,t)} are correlated Wiener processes. Their increments are 
expressed as 

M / M \ M / M \ 
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where {Vi(i)} and {W^ a \t)} are independent Wiener processes. The statistics of {d^ a \<j>, i)} are specified by 
(dC (a) (0,t)) = and (dC (Q) (0,t)dC (/3) (0^)> = C {a ^{4>)dt, where C^\<j>) is a # WxAr correlation matrix defined as 



M / M \ I M 



fc=l \i=l / \i=l 

M / M \ I M 

+ e E E^ (a) ) ff ^ (a) ) E^^ (/3) )^^ (/3) ) ] <w- ( 4 ) 

fe=l \i=l J \j=l 

Note that C^ a, ^\<p) 1S periodic in for all a, and its (a, /3)-component depends only on (j)^ and Since 
C( a '^(4>) is a symmetric positive semi-definite matrix, we can also express d^ a ^(4>,t) using N independent Wiener 
processes {WW(t)} as dC, {a \4>,t) =Y J N p=1 B ( - a - p \<t>)dW^\t), where is a real symmetric matrix satisfying 

Ef=i B^ a ^(4>)B^^(cf)) = C^\(t>). The transformed drift coefficients can be calculated as 

A(a) W=» + lo^> C(a>a) (cl>), (5) 

where we utilized the fact that the right-hand side of Eq. ([2]) depends only on in calculating the Wong-Zakai 
correction term [To| . The original N vector Stratonovich Langevin equations ([T]) with N + 1 independent vector noises 
and {r7 (Q) (i)} are now reduced to AT scalar Ito stochastic differential equations with N correlated scalar noises 
{d(( a ^ (4>, t)}. The corresponding Fokker-Planck equation (FPE) describing the evolution of the probability density 
function (PDF) P(4>, t) of the phase variables is given by [10( 

a N a 1 N N n2 

a=l r a=l/J=l r ^ 

We now invoke the averaging approximation Q to this FPE. We introduce new slow phase variables i/> = 
• • • , V (Ar) ) as <j)^ =ut + tfj^ (a = 1, ■ • ■ , iV), and their PDF 

Q(i/>, i) = Q({^ Q >}, t) = P({^ Q > = art + </> (a) }, t). (7) 

With sufficiently weak external noises, Q varies slowly compared with the oscillator natural period, T = 2tt/lu. We 
can thus average the drift coefficients A^ a \cf)) and the diffusion coefficients C^ a, ^(cf)) of the FPE over the period T 
keeping Q constant. The resulting averaged FPE for Q is given by 

N N 

1 1- 

(8) 



a=l (3=1 Y r 



The drift coefficient A^((p) simply yields lj after averaging due to the periodicity of C a ' ,9 (0) in <j)( a \ which vanishes 
in the new variables. The averaged diffusion coefficients D^ a ^\xjj) are given by 



1 f t+T 

£)(«*.« (^) = - J C (a ^({^ = cot'+^})dt' = Dg(^-^) + eh{0)S a ,p, (9) 

where we utilized the fact that C^ a, ^\4 > ) depends only on cf>^ and <f>^\ and introduced the correlation function g(0) 
of Zi(4>) and Gik(4>) as 

and similarly the correlation function h(9) of Zi{4>) and Hik(4>) as 

h ( ) = E H^Hiki^Zj^ + 6)H jk {<t>' + 0)d<t>'. (ff) 

27r J -«i,j,k=i 
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Clearly, g(0) > and h(0) > (we exclude the non-physical case Z(<p) = const.). Using the periodicity of Z(<f>) 
and G(<p) in <p, it can also be proven that g(9) = g{—9) and g(0) > g{9). Since Z{<f>) and G(<p) are smooth functions 
of (j), g(9) has a quadratic peak at 9 = 0. It can also have other quadratic peaks at 9 ^ 0, e.g. 9 = ±27r/3, when .£(</>) 
contains non-negligible high-order harmonics or when the common noise is introduced multiplicativcly. 

To analyze the phase relationship between the oscillators, we focus on the PDF of the phase difference. Without 
loss of generality, we first introduce the two-body PDF of and V> (2) as ^ (2) ,t) = J <# (3) • • • # (Ar) <2(V>,*)- 

The evolution equation for R(ip( x \ ip( 2 \t) can be derived from Eq.© by integrating over all other phase variables as 

!*»»,*«,.) = I W ,(O) + e M0)){(^) 2 + (^) ! j fl+s? ^(l, 9W <..-^.,«). (12, 

Furthermore, by transforming the two phase variables to the mean phase and the phase difference, ip = (t/n 1 ) + tp^)/2, 
9 — — tp( 2 \ the above equation can be further decoupled as 



J^foM) = \{D[g(0)+ g (9)]+eh(0)}-^S(^t), 

Q-U(6,t) = — {D[g(0)-g(9)]+eh(0)}U(9,t), (13) 

where S (ip , t)U (9 , t) = R(ip^ = tp + 9/2, -p^ = ip - 9/2, t). It is clear that Eq. JT3D has a unique final stationary 
state, where the PDF of the mean phase ip is uniform over the limit cycle, So(ip) = l/2ir, and the PDF of the phase 
difference 9 is given by 

W) = D{g(0)-g(9)] + eh(0y (M) 

where uq is a normalization constant. 

We now examine the consequences of the above results. Our argument holds generally for arbitrary g(9) that 
satisfies our assumptions, namely, for a general class of limit-cycle oscillators. When only the independent noises are 
given, D = and e > 0, Uo(9) is simply uniform, so that the oscillators are completely desynchronized. When only 
the common noise is given, D > and e = 0, Uq(9) diverges at 9 = while remaining positive because g(0) > g(9), 
so that the phase difference between any pair of oscillators accumulates at zero, resulting in noise-induced complete 
synchronization. As e is increased from zero, Uq(9) becomes broader, but its peak at 9 = remains as long as D > 0, 
i.e. the oscillators still concentrate coherently around 9 = 0. As we mentioned previously, g{9) may have multiple 
peaks in addition to 6 = 0. Then, the above discussion also holds for such values of 9. Multiple peaks of g{9) lead 
to the clustering behavior of the oscillators, a well-known phenomenon in coupled oscillators [9(, but in the present 
case, it is caused by the combined effect of the phase sensitivity and the common noise alone. More generally, Uq(9) 
can exhibit a wide variety of non-uniform "coherent" distributions depending on the functional form of g{9). 

We can also examine the statistical stability of the synchronized state 9 = and the dynamics of 9 
around it. From Eq. (fT3"]) . we obtain the corresponding Ito stochastic differential equation for 9 as d9(t) = 

V2{D[g(0) - g(9)] + eh(0)}* dw(t), where w(t) is a Wiener process. Focusing on the region around 9 = 0, we 
approximate g(9) around its 9 = peak as g{9) ~ g(0) — (l/2)\g"(0)\9 2 , utilizing the facts that g'(0) = and 
g"{0) < 0, where ' denotes d/d9. We then obtain 



d9(t) ~ ^D\g"(0)\9(t)d Wl {t) + y/2eh(0)dw 2 (t), (15) 

where the noise is decomposed into multiplicative and additive parts using two independent Wiener processes wi^{t)- 
This is simply a linear random multiplicative process with an additive noise pT| . Let us ignore the additive noise 
diU2(t) for the moment. Using the Ito formula [10(, the equation for the logarithm of the absolute phase difference is 
obtained as dln\9(t)\ = -±D\g"(0)\dt + y /D\g"{0)\dw 1 (t), so that the average Lyapunov exponent of the completely 
synchronized state 6 = is given by A = — ^D\g"(0)\ < 0, which is always negative, i.e. 9 — is always statistically 
stable. When the common noise is additive, G(cp) is a constant matrix, and we recover the previous results 0, lif- 
When weak independent noises exist, 10 mostly remains small but occasionally exhibits large bursts, a typical behavior 
known as noisy on-off intermittency [11 1. We then expect a power-law PDF of the inter-burst intervals of 9(t) with an 
exponent —1.5, and also a power-law PDF of the phase differences around 9 = 0, whose exponent is always —2 in the 
present case [ll| (results not shown; see Ref. 0|). When g{9) has multiple peaks, we can estimate the stability and 
fluctuations around the other peaks in a similar fashion, and we expect intermittent transitions between the clustered 
states 1141 . 
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We now demonstrate the noise- induced synchronization and clustering numerically. As the first example, we consider 
uncoupled Stuart-Landau (SL) oscillators, X = (x,y), F(X) = ( x — c a y — (x 2 + y 2 )(x — c 2 y), y + c x— (x 2 + y 2 ){y + 
C2x) ), subject to independent additive noises, H{X) = diag(l, 1), and to the following four types of additive or 
multiplicative common noises, G\(X) = diag(l,l), G 2 (X) = diag(x, y), Gs(X) = diag(l + 4xy,0), and G±(X) = 
diag(x, xy). The SL oscillator is the simplest limit-cycle oscillator derived as a normal form of the supercritical 
Hopf bifurcation We fix the parameters at Co = 2 and C2 = — 1, with which the natural frequency becomes 
lu = Co — C2 = 3. The phase sensitivity function is analytically given as Z{<f>) = V2( sin (cf> + 3vr/4) , sin (0 + tt/4) ) 0]. 
From Eq. pTj) . we obtain the corresponding correlation functions as gi(9) = 2cos6', 32 (^) = cos 2 0, .93(6*) = cos3#, 
g±{0) = (cos# + 8cos 2 + cos30)/16, and h(0) = 2, from which we can calculate Uq(0). We thus expect noisy 
synchronization (1-cluster), 2-cluster, 3-cluster, and intermixed coherent distributions of 9 to be observed. Figure [T] 
compares the results of direct numerical simulations using N — 200 oscillators with the analytical results, where the 
noise intensities are fixed at D = 0.002 and e = 0.0001. To realize the Stratonovich situation, the numerical simulations 
are performed using colored Gaussian white noises generated by the Ornstein-Uhlenbeck process rz(t) = —z(t) 
with a small correlation time r = 0.05, where £(i) is a Gaussian white noise of unit intensity [Icj . As expected, various 
synchronized or clustered states are realized, and their PDFs are fitted nicely by the theoretical curves. 

As the second example, we consider uncoupled FitzHugh-Nagumo (FN) oscillators 8], X = (u,v), F(X) = 
( e(v + c — du), v — v 3 /3 — u + I), subject to either an additive or multiplicative common noise, G±(X) = diag(0, 1) 
or G2(X) = diag(0, v), and also to an additive independent noise, H(X) = diag(0, 1). The noises are applied only to 
the variable v corresponding to the membrane potential. Fixing the parameter values at e = 0.08, c = 0.7, d = 0.8, 
and / = 0.875, the limit cycle becomes symmetric with a natural frequency of u> ~ 0.1725. The phase sensitivity 
function Z v (<f>) can be calculated numerically using the method devised in 0, @]. Figure [5] compares the results of 
direct numerical simulations with the analytical results at D = 0.005 and e = 0.0005 using N = 200 oscillators. Either 
synchronized or 2-cluster states are realized for the additive or multiplicative common noises, and their PDFs are well 
fitted by the theoretical curves calculated using the numerical Z v (<fi). 

Summarizing, we developed a global formulation of synchronization and clustering phenomena in ensembles of 
uncoupled limit-cycle oscillators induced by a common noise. The common noise acts as a state-dependent noise 
on the phase difference, which yields the ^-dependent effective diffusion constant for U(8,i) in Eq. (TT3"1) , an d results 
in the non-uniform coherent stationary distribution Uq(0) in Eq. (|14[) . In our formulation, the combination of the 
common and independent noises is a special case of more general correlated noises, and the synchronized or clustered 
state is a special case of non-uniform coherent distributions. Thus, we can generalize the notion of common-noise- 
induced synchronization to correlated-noise-induced coherence. This insight will be helpful in understanding various 
noise-induced synchronization phenomena. 
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FIG. 1: Stuart-Landau oscillators, (a) Synchronized, (b) 2-cluster, (c) 3-cluster, and (d) intermixed states. The insets display 
instantaneous distributions of the oscillators on the limit cycle. 




FIG. 2: FitzHugh-Nagumo oscillators, (a) Variable v(<f>) and phase sensitivity function Z v (<f>), (b) correlation functions 51,2(6') 
calculated from Gi.2(X) and Z v ((j>), (c) synchronized state, and (d) 2-cluster state. The insets display snapshots of the 
oscillators. 



